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Abstract: Nulling interferometry is an astronomical technique that 
combines equal wavefronts to achieve a deep rejection ratio of an on-axis 
star, and that could permit to detect Earth-like planets in the mid-infrared 
band [5 - 20 jum]. Similarly to what is done in the near-infrared, high 
frequencies spatial filtering of the incoming beams can be achieved using 
single-mode waveguides operating in the mid-infrared. An appreciable 
reduction of the instrumental complexity is also possible using integrated 
optics (IO) devices in this spectral range. The relative lack of single-mode 
guided optics in the mid-infrared has motivated the present technological 
study to demonstrate the feasibility of dielectric waveguides functioning at 
longer wavelengths. We propose to use selenide and telluride components 
to pursue the development of more complex IO functions. 
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1. Introduction 

In the last ten years, the photonics field has provided a number of valuable technological solu- 
tions to astronomers to improve the performance of the observatories instruments. In particular, 
optical fibers and integrated optics (IO) devices have already been part [1, 2], or will be part 
in the future [3, 4] of interferometric facilities. As a result, this context has raised the question 
of extending the use of photonics devices to the Darwin Infrared Space Interferometer of the 
European Space Agency (ESA) [5]. 

This mission, based on the implementation of nulling interferometry [6], is devoted to the 
search and characterization of terrestrial exoplanets in the mid-infrared range thanks to its com- 
bined high contrast and high spatial resolution (long -baseline interferometry) features [7]. The 
principle of a two aperture nulling interferometer is to recombine artificially n phase-shifted 
wavefronts from an on-axis star in order to achieve their mutual cancellation. The signal from 
an off-axis object - like a planet orbiting the star - is transmitted through the instrument by tun- 
ing the interferometer baseline so that a total phase shift of 2 n is achieved. A key requirement of 
nulling interferometry is equal wavefronts. Reducing the wavefront errors (e.g. pointing errors, 
optical defects due to polishing and coatings) can be obtained with single-mode waveguides 
as modal filters [8]. Using single-mode Integrated Optics (IO) components merges the func- 
tions of modal filtering and beam combination, providing a stable and compact instrument with 
regards to external mechanical and thermal constraints. To date, the implementation of modal 
filtering with single-mode waveguides has become the major option for space-based nulling 
interferometry. 

The lack of single-mode waveguides for the mid-infrared in the photonics field has pushed 



advanced technological studies to develop such devices [9, 10, 11]. In this context, we have 
initiated, under ESA contract, the IODA program (Integrated Optics for ZMrwin) to investigate 
the feasibility of a single-mode Integrated Optics concept operating over the full 5-20 jxm 
spectral range. Our program has focused on the feasibility of an IO beam combiner based on 
either conductive waveguides [12] or on dielectric materials. 

In this paper, we present results obtained using selenide and telluride materials to manufac- 
ture planar waveguides for the mid-infrared range, the first step toward further efforts for the 
development of channel waveguides. We show that the modal behavior of the various structures 
can be experimentally investigated and assessed with the m-lines method, known in the visible 
and near-infrared range, and which we have successfully transferred to the mid-infrared range. 
To our knowledge, this paper presents the first study of the modal properties of selenide and 
telluride thick films at longer wavelengths. 

2. Investigating the properties of slab waveguides through m-lines study 

2.1. Short theoretical analysis on fundamental parameters 

It is well known that a planar waveguide must fulfill the total reflection condition at the interface 
between film-substrate and film-superstrate [13], corresponding to « c >max(n su b, n S up) with 
n c , «sub and « sup the refractive indices of the film, substrate and superstrate, respectively. To 
propagate at least one mode, the thickness d of an asymmetric planar waveguide must be greater 
than the cut-off thickness of the fundamental mode do given, in a more general form, by: 



where k=2n/X is the wave vector, m the mode order starting from 0, and g a polarization 
factor with g=l for TE polarization and g=n\jn\ vh for TM polarization. Thus, the quantity d m 
increases with the operating wavelength and also depends on the dispersion of the material. 
Here, we briefly evaluate from a simple theoretical approach the possibility to benefit of a 
single-mode planar waveguide over the full Darwin range (5 - 20 jJ.m) using As2Se3 glass 
as the layer and AS2S3 glass as the substrate. The superstrate is air. We assume that these 
materials are transmissive over the entire range of interest, neglecting possible absorption lines. 
We also consider only the case of TE polarization with g=l. We derive from the refractive 
indices of these bulk materials [19] a simple fit of the dispersion laws which are given, to a 
first approximation, by n c (X) = -0.00605 xA+2.85091 and « sub (A)=-0.0095 xA+2.47128 for 
As2Se3 and AS2S3 respectively, with X given in microns. The adopted models are justified 
by the fact that the refractive indices of these materials decrease almost linearly in the 5-20 
jj.m range. Under those assumptions, we compute the cut-off thicknesses do and di for modes 
m=0 and m=l and for 5 /im<A<20 fim. The acceptable slab thickness d for single-mode 
behavior over our spectral band must comply with c/o(20/zm)<c/<rfi(5/im). However, since 
we found Jo(20/im)=2.903 fim and iii(5/im)=2.577 jim, the previous condition cannot be 
fulfilled. The wavelength X\ for which do{X\)<d\{5pLn\) is found to be X\=\l .2 fim. Thus, it 
theoretically exists a film thickness, e.g. d=2.54 jj,m, for which a planar waveguide following 
these dispersion laws propagates the fundamental mode only over a 5 - 17 jj,m band, i.e. 
almost the full expected band. 

We underline that this modal analysis is based only on the study of the opto-geometrical 
parameters of the waveguide and is thus sensitive to the numerical fit. To better constrain 
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Fig. 1. Typical prism coupler setup. When we are in the condition of total reflection at 
the base of the prism (0 p smaller than critical angle 9 C ), the modes of the slab waveguide 
can be excited through the evanescent field existing in the air gap. The air gap thickness is 
controlled with a pushing screw at the back of the waveguide. A and B are the angles of the 
prism. 6i m and ftj.m are the algebraic angles with respect to the normal of input and output 
prism facets corresponding to the excitation of mode m. 8 m is the angle of the propagating 
ray corresponding to mode m. 

the trade-off between fundamental mode confinement and higher-mode filtering properties 
(which requires to be far from the cut-off), a more detailed electromagnetic approach through 
Beam Propagation Method (BPM) simulations will be required. However, our results show 
that a good compromise could reasonably be reached by splitting the 5-20 jUm band into no 
more than two sub-bands, with two dedicated waveguide designs, and on the condition that a 
sufficient index difference between the film and the substrate is ensured. 

2.2. Thick film deposition 

Extending the 10 solution to the mid-IR requires infrared materials with suitable features in 
this spectral domain. Chalcogenide glasses are interesting candidates, due to their good trans- 
parency in the mid-infrared, their high refractive index, and the relatively ease of preparation 
in films form. We have selected As2Se3, Te2As3Se5 and TeAs4Ses chalcogenide compositions 
to undergo thick film deposition on an AS2S3 substrate. Glasses containing elements like sul- 
phur (S), selenium (Se) and tellurium (Te) transmit at longer wavelengths: sulphide glasses (i.e. 
AS2S3) transmit up to 12 jim, selenide glasses (i.e. As2Se3) transmit up to 16 /im, and tel- 
luride glasses transmit up to 20 /im [14], with an average transmission of about 35% for the Te 
glass. Note that this transmission value includes the Fresnel losses at the interfaces due to the 
high refractive index of the glass. The deposition process used is thermal evaporation, which 
reduces the column-like effects in the structure and provides dense films with higher refractive 
index with respect to RF-sputtering techniques. The selenide films were deposited from As2Se3 
commercial bulk glass, while the telluride films used bulk glass that was synthesized in the 
laboratory from commercial components like Tellurium (Te), Arsenic (As) and Selenide (Se). 

The strategy adopted to develop the chalcogenide planar waveguides is the following. For 
both selenide and telluride compositions, we have selected AS2S3 chalcogenide commercial 
glass as the substrate - the superstate is air -, which presents the lowest refractive index 
among the different glasses available («=2.46 at A=1.2 /im and «=2.38 at A=10 /im). As under- 
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Fig. 2. Experimental arrangement of the mid infrared m-lines experiment: QW, quarter- 
wave plate; M, mirror; BS, beam-splitter; A, absorbing refractory brick; P, polarizer; B, 
beam-expander; R, reticle; F, focusing optics; S, system prism - film; C, infrared camera. 



lined in Subsection 2.1, the challenge consisted in obtaining the appropriate index difference 
An = n c — n su b and thickness d in order to propagate at least one mode. Eq. (1) shows that, for 
a given An, d m is proportional to X. Thus, at longer wavelengths the target thickness must be 
significantly higher than the usual thicknesses (« 0.5 jJ,m) deposited for near-IR applications. 
Because of the novelty of the manufactured samples, an unknown is put on the achievable index 
difference and thickness. Consequently, we have experimentally investigated the modal behav- 
ior of the sample, the refractive index and the thickness of the films through the mid-infrared 
m-lines experiment developed at LAOG (Laboratoire d' Astrophysique de Grenoble). 

2.3. Experimental technique 

The m-lines method, based on prism coupling theory, was first proposed by Tien et al. [15] 
and is commonly used as a characterization method for telecommunication components in the 
visible and near-infrared range. In this study, we have implemented this method at A=10.6 fim. 
In this technique, a prism is brought in close contact to the surface of a planar waveguide. 
Using a pushing screw placed under the waveguide (see Fig. 1), a small and localized optical 
contact is created between the prism base and the film. This provides a thin air gap that should 
optimally be a quarter of the operating wavelength to maximize the coupling effect. When a ray 
reaches the base of the prism with an incident angle smaller than the critical angle defined by 
the prism index, it becomes totally reflected. However, the overlap of the evanescent fields in 
the air gap as predicted by electromagnetic theory [16] implies that light can be coupled into the 
waveguide similarly to a tunneling effect. The existing modes of the waveguide will be excited 
if the phase matching conditions given by Eq. (2) are fulfilled 



n p cos(0 p ) = n c cos(0 m ) (2) 

where n p and n c are respectively the prism and the thick film indices, and p and 9 m are the 
angular directions of the rays as illustrated in Fig. 1 . Exciting the low-order modes of the struc- 
ture requires a prism with higher index than the slab waveguide index. Under these conditions, 
and for the incident rays that fulfill the phase conditions, the coupling phenomenon will lead to 
the observation of "black lines" at discrete positions 02. m m the reflected beam pattern. These 




Fig. 3. Output of the m-lines experiment in the mid-IR (left) and in the near-IR (right). The 
image of the reticule is used to point the black lines by turning the prism stage. The mid-IR 
pattern presents a circular diffraction pattern due to the high coherence of the laser source. 

lines are characteristic of the different modes of the slab waveguide. From the measurement of 
the angles 02. m, the effective index « e ff,m of mode m is computed through: 



where B is one of the prism angles. The computation of the thick film index « c and its thickness 
d is obtained by implementing a non-linear least squares method (see Subsection 2.4). 

The experimental setup to measure the mode indices at 10.6 flm is shown in Fig. 2. The 
infrared source is a CO2 laser. To probe either TE or TM modes, the quarter-wave plate QW 
modifies the intrinsic linear polarization of the source into an elliptical one, then the grid po- 
larizer P selects the desired polarization. The Zinc Selenide (ZnSe) optics B and F produce a 
converging beam which is focused on the optical contact. This beam includes the angular direc- 
tions that fulfill Eq. (2), simplifying the search for mode lines. The reticule is used to align the 
normal of the input face of the prism with the optical axis of the bench in order to solve Eq. (3). 
The black lines in the reflected beam are observed with a mid-infrared camera (see Fig. 3). 

For the experiment at 10.6 jJ.m, we used a high-index germanium prism with an angle 
B=30°±0.01° and a refractive index « p =4.0±0.001. In one case, the m-lines experiment was 
run at shorter wavelengths (A=1.196 jim) and for which we used a silicon prism with an angle 
B=44.8416°±5x 10 -4 and a refractive index « p =3.5193±10 -4 . In our system, the position of 
the black line is measured with an accuracy of ±0.01°. 

2.4. Data processing 

In this paper, we have assumed that the tested samples are asymmetric slab waveguides (i.e. 
Air/As?Se3/As?S3) with a step index distribution. To obtain a correct numbering of the modes, 
we used the property of linearity between the squared spectrum of mode indices n% s (m) and 
the mode numbers set (m+1) 2 , as formalized by Lee et aZ.[17]. Thus, when more than three 
experimental mode indices are available, their numbering is derived from the computation of 
the linear regression coefficient p of the experimental curve [« 2 ff (f7z),(m+l) 2 ] and compared to 
unity (see Tables 1,2). In this paper, we assume that mode numbering starts with 777=0. Succes- 
sively, we compute n c and d by numerically solving [18] the waveguide dispersion equation for 
a step index distribution planar waveguide given by Eq. (4) [13]: 
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a 


N/A 


2.7719 


N/A 
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+25.65° 
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2.7697 


2.7702 


0.0005 
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+25.50° 


+0.15 


2.7682 


2.7679 


0.0003 
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+25.21° 


+0.29 


2.7654 


2.7649 


0.0005 
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+24.81° 


+0.4 


2.7614 


2.7613 


0.0001 
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+24.39° 


+0.42 


2.7573 


2.7570 


0.0003 
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+23.79° 


+0.6 


2.7513 


2.7519 


0.0006 
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+23.25° 


+0.54 


2.7459 


2.7462 


0.0003 


9 


+22.67° 


+0.58 


2.7400 


2.7399 


0.0001 


Refractive index n c 2.7732 ±1.3 x 1(T 4 


Thickness d 13.74 ± 0.05 jlm 



Table 1. TE mode indices obtained experimentally on sample As-1 at A=1.196 flm with 
the m-lines setup. The error on the measured angle is ±0.01°. The error on experimental 
n e ff is ±1.3x 10~ 4 . We found a maximum coefficient p=-0.99923 for the presented mode 
numbering. The silicon prism has n p =3.5193 and B=44.8416°. The substrate refractive 
index at A=1.196 jlm is 2.47[19]. Note that the two first modes could not be experimentally 
observed with the silicon prism and are labelled as "a" for "absent". The sign "N/A" means 
"not applicable". 
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To better constrain the result given by the linear regression coefficient method, we use the 
derived quantities n c and d to compute the theoretical mode indices for different mode num- 
bering by analytically solving Eq. (4). The correct mode numbering is then obtained through a 
best fit between the theoretical and experimental distributions of indices. 

3. Results and discussion 

3.1. Selenide components 

The samples, with a target thickness around 10 jttm, were characterized in the near-IR (A=1.2 
|im) and in the mid-IR (A=10.6 jlm) with the m-lines setup presented in Subsection 2.3 . Two 
samples were produced and respectively referenced as As-1 and As-2. Sample As-1 was char- 
acterized both at A=1.2 fim and at A=10.6 fim, while sample As-2 was characterized in the 
mid-IR only, Including the near-IR characterization is advantageous, since a multimode behav- 
ior is expected from the moment the film thickness is greater than 1 jlm. Thus, this increases 
the number of mode lines that can be detected and measured. 

Table 1 reports the TE mode indices and the derived parameters for As-1 at X= 1.196 |im. 
Sixteen mode lines were experimentally observed with the sample As-1, but only the first eight 
modes are presented here. Note that measurements of TM modes were not carried out because 
the two optical components QW and P in Fig. 2 were not available at the time of the corre- 
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+55.04° 


N/A 


2.667 


2.6666 


0.0004 


1 


+46.89° 


8.15 


2.599 


2.5992 


0.0002 


2 


+36.09° 


10.8 


2.488 


2.4868 


0.0012 


Refractive index « c 


2.689 ± 0.001 






Thickness d 


13.27 ±0.03 fim 






sample As-2 





+48.58° 


N/A 


2.614 


2.6138 


0.0002 


1 


+42.61° 


5.97 


2.557 


2.556 


0.001 


2 


+33.61° 


9.00 


2.4601 


2.4606 


0.0005 


Refractive index n c 


2.633 ±0.001 






Thickness d 


14.49 ± 0.03 jum 







Table 2. TE mode indices obtained experimentally on samples As-1 and As-2 at A=10.6 
jlm. The error on measured angle is ±0.01°. The accuracies on prism index and angle lead 
to an error on experimental n e ff of ±10~ 3 . The maximum coefficient is p=-0.99982 for the 
given mode numbering. The prism has n p =4.00±0.001 and B=30.0±0.01°. The substrate 
refractive index at A=10.6 flm is 2.38. The sign "N/A" means "not applicable". 



sponding tests. Thus, only the intrinsic linear polarization of the source, corresponding to TE 
modes, was used for those two samples. The derived quantities are n c =2. 7732 ± 1.3 xlO" 4 and 
of=13.74 ± 0.05 jiim. A similar characterization was made at A=1.55 flm, in order to observe 
a chromatic change in refractive index. This led to values « c (1.55/zm)=2.7407±1.3x 10^ 4 and 
d=13. 16 ± 0.03 jim. We observe a change in refractive index of approximately 0.04 as well as a 
small variation of the film thickness. Note that the m-lines experiment provides a measurement 
of the pair of quantities (n c ,d) for the film area where the optical contact is obtained. When 
shifting from k=\ .2 jj,m to A=1.55 /im, we also moved the prism on the film, which has an 
effect on the measured film thickness. 

We performed M-lines characterization in the mid-infrared at A=10.6 |im for the two sam- 
ples As-1 and As-2. At longer wavelengths, the germanium prism described in Subsection 2.3 
was used. Table 2 shows the experimental results obtained with the two samples. The value of 
the substrate index, « s =2.380, has been taken from the literature [19]. 

For the two samples, three mode lines have been observed with our setup. The derived theo- 
retical mode indices are very close to the experimental ones. The dispersion of the mode indices 
is low, which implies a good homogeneity of the index distribution. The derived refractive in- 
dices are « c =2.689 ± 0.001 for As-1 and « c =2.633 ± 0.001 for As-2, which shows that the 
densification of the layer was lower during the second run. The sample thickness is approxi- 
mately 13 jJ.m for As-1 and is consistent with the value obtained at shorter wavelengths. The 
thickness is slightly higher for As-2 with a value of approximately 14.5 fxm. 

We measure mid-IR refractive indices significantly lower than in the near-IR. This is consis- 
tent with the dispersion law of the material, which predicts a decrease in the refractive index 
with increasing wavelength. We also note that the achieved thicknesses are greater than the 
targeted ones, which shows the deposition process needs to be better stabilized. Finally, a mul- 
timode behavior of the planar structures is clearly identified in the near and mid-IR ranges. 
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0.006 


2 


+28.13° 


+15.43 


2.394 


2.395 


0.001 
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a 


N/A 
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2.686 


N/A 
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+45.26° 


N/A 


2.583 


2.5887 


0.0057 
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+30.13° 


+15.13 


2.419 


2.4298 


0.0108 


Refractive index n c 




2.718 ±0.001 




Thickness d 




10.69 ± 0.03 Aim 





Table 3. TM (up) and TE (down) mode indices obtained experimentally on sample TAS-1 
at A=10.6 flm. The error on measured angle is ±0.01°. The accuracies on prism index and 
angle lead to an error on experimental n e ff of ±10~ 3 . The same germanium prism as in 
Table 2 is used. The sign "a" means "absent" and "N/A" means "not applicable". 

3.2. Telluride components 

In this approach, the thick films were deposited by evaporation on a AS2S3 substrate from 
Te?As3Se5 and TeAs4Ses bulk glasses made by LPMC (Laboratoire de Physico-chimie 
de la Matiere Condensee). Two samples referenced TAS-1 (TeAs4Ses/As2S3) and TAS-2 
(Te?As3Se5/As2S3) were produced by thermal evaporation and characterized for TE and TM 
polarizations at A=10.6 jim. The same germanium prism was used to achieve light coupling. 
The experimental measurements on TAS-1 are reported in Table 3. For TAS-2, we present the 
derived refractive index and thickness only in Table 4. Mode lines have been observed for the 
two samples. The numerical analysis confirmed that the dark lines correspond to propagating 
modes. In this case as well, more than one line were observed. 

For sample TAS-1, the fundamental mode in TE polarization could not be observed, oppo- 
sitely to the TM case. For each mode, the measurement of the line angular position for the 
two polarizations was done successively. In this way, we could experimentally observe that 
there was a slight angular shift for the same mode between the two polarization. The ra-lines 
measurement was done for TE and TM polarizations at the same point of the film. The results 
show a small difference between TM and TE index n c , which supposes an hypothesis of bire- 
fringence of the film. For this first telluride component, an approximate thickness of the film of 
about 10 /xm has been achieved. 

The sample TAS-2 was characterized with a similar procedure. Three mode lines have been 
experimentally observed but only in TE polarization. The corresponding TM case could not 
be observed despite our efforts in optimizing the coupling efficiency. For TAS-2, we therefore 
present the derived refractive index and thickness only in Table 4. 



sample 


TAS-2 (TE polarization) 








Mode 


Measured 


Angle 


Observed 


Theoretical 


Diff . 




angle 


separation 


nrff 


n e ff 







+73.59° 


N/A 


2.772 


2.776 


0.004 


1 


+52.45° 


+21.14 


2.647 


2.640 


0.007 


2 


+30.00° 


+22.45 


2.417 


2.419 


0.002 


Refractive index n c 




2.821 ±0.001 




Thickness d 




8.79 ± 0.03 jum 





Table 4. TE mode indices obtained on sample TAS-2 at A=10.6 flm. The error on measured 
angle is ±0.01° and the accuracies on prism index and angle lead to an error on experimen- 
tal n eff of ±10~ 3 . 

Compared to sample TAS-1, TAS-2 presents a higher refractive index at A=10 flm. The 
possibility to obtain a high refractive index is an important result if a "full-telluride" (film and 
substrate) is foreseen. Using a telluride glass as a substrate will indeed raise the question of 
the achievable index difference between the film and the substrate, since telluride bulks are 
expected to present higher index compared to selenide bulks. 

3.3. Impact of the experimental results on the requirements for single-mode slab waveguides 

These experimental results show that the refractive index obtained with As2Se3 under thick 
film form is lower (by approximately 0.15) than in bulk form [19]. Since the index differ- 
ence between the film and the substrate is an important parameter, we have now a less fa- 
vorable case compared to the analysis of Subsection 2.1. Deriving the new model n c (X) 
= -0.0127 xA+2. 78532 from the data of Table 1 and 2, we find that we would need two 
As2Se3/As2S3 slab waveguides to cover the 5-20 flm band: the first one with a thickness 
d=3.\l flm for the sub-band 5 - 12.5 |im, and the second one with a thickness d=7.0 flm for 
the sub-band 12.5 -20 |im. 

Extending the first single-mode sub-band can be achieved with the telluride glass Te2As3Ses 
as the guiding layer, since the experimental results of Table 4 show that the achievable refrac- 
tive index is higher than for As2Se3. 

However, the AS2S3 substrate transparency remains an issue - it does not extend up to 20 
flm. Therefore, covering the whole 5-20 flm band would require an examination of a 20 flm 
transparent telluride substrate, whose refractive index is sufficiently low to ensure a An of about 
0.42. 

4. Conclusions 

This paper has discussed the optical characterization of chalcogenide-glass-based thick films 
manufactured for the purpose of future mid-IR integrated optics development. We have fully 
developed an ra-lines setup operating at A=10.6 flm based on prism coupling theory. This setup 
was used to characterize the modal behavior of the first mid-infrared planar waveguides by ob- 
serving mode lines. In this study, two different compositions were successfully tested at 10 jiim 
with the first observation of three mode lines for TE and TM polarizations. A simple numerical 
analysis based on a best fit between the experimental and theoretical mode indices distributions 
has permitted us to derive the refractive index and thickness of the different samples. Our results 
have underlined the possibility of producing full-selenide multimode planar waveguides with a 
sufficient film thickness. The near and mid-IR results are consistent, since similar thicknesses 
are obtained in the two spectral ranges. Extending the transmission range of the components to 



20 jUm may be possible with the telluride solution. 

Following these results, a complementary technological study has been conducted to in- 
vestigate the feasibility of a first full-telluride single-mode planar waveguide. This study has 
produced positive results. As a consequence, the question of etching of a channel waveguide to 
be used as a modal filter for nulling interferometry has been investigated straightforward and 
the results recently obtained will be part of a forthcoming paper. 
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